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ABSTRACT 

Several recent measurements have been made of the angular correlation function 
w(9, m) of faint galaxies in deep surveys (e.g., in the Hubble Deep Field, HDF). Are the 
measured correlations indicative of gravitational growth of primordial perturbations or 
of the relationship between galaxies and (dark matter dominated) galaxy haloes? A first 
step in answering this question is to determine the typical spatial separations of galaxies 
whose spatial correlations, £(r, z), contribute most of the angular correlation. 

The median spatial separation of galaxy pairs contributing to a fraction p of the 
angular correlation signal in a galaxy survey is denoted by (SQ) and compared with 
the perpendicular distance, rj_, at the median redshift, -2 me d, of the galaxies. Over a 
wide range in spatial correlation growth rates e and median redshifts, is no more 
than about twice the value of r±, while is typically four times r±. 

Values of for redshift distributions representative of recent surveys indicate that 
many angular correlation measurements correspond to spatial correlations at comoving 
length scales well below l/i -1 Mpc. For = 1 and Ao = 0, the correlation signal at 4" 
predominant in the [Villumsen et al.| (1996) estimates of w(9, m) for faint HDF galaxies 
corresponds to (HDF) « 40/i _1 kpc; other cosmologies and angles up to 10" can 
increase this to (HDF) <,200/i _1 kpc. The proper separations are (1 + z) times 
smaller, where 1 <, z <, 2. 

These scales are small: the faint galaxy angular correlation measurements are at 
scales where halo and/or galaxy existence, let alone interactions, may modify the spatial 
correlation function. These measurements could be used to probe the radial extent of 
haloes at high redshift. 



Subject headings: cosmology: theory — galaxies: formation — galaxies: clusters: general — 
galaxies: distribution — cosmology: observations 



1. INTRODUCTION 

The growth of primordial density fluctuations into 
galaxies and clusters of galaxies is an essential ele- 
ment of understanding how structure forms in the 
Universe. A common way to statistically represent 
this structure from observed galaxies is to calculate 
the two-point auto-correlation function, £(r, z), which 
can be approximately parametrised as a power law (in 
spatial separation of galaxy pairs) which increases in 
amplitude as a function of time, according to a single 
parameter e, (e.g., Groth fc Peebles] 1977): 

(1) 



f(r,z) = (r /r)^(l + z)-( 3 + e -^ 



where r and rg are expressed in comoving coordi- 
nates and 7 represents the approach to homogeneity 
at larger length scales. 

Observed values are typically ro ~ 5/i _1 Mpc and 
7 « 1-7 — 1-8 (e.g., Davis & Peebles 1983; Loveday 



t al. 1992) for the low rcdshift general galaxy pop- 



ulation. 



Davis . 



Peebles' (1983) analysis is consistent 

1 Mpc, 



with this power law on scales lOh 1 kpc , 



10h- 



though the correlations appear somewhat stronger at 
the 



The validity or otherwise of the high observational 
estimates of e or of major changes in the general 
galaxy population is not the subject of this article. 
In fact, as will be defined and seen below, the param- 
eter of interest r p ^ is only weakly sensitive to e. 

The motivation for this article stems from the at- 
tempts to indirectly measure the dependence of £ on 
z. Since estimation of spectroscopic redshifts requires 
many more photons than the photometric detection of 
a galaxy, the projection of £ onto the celestial sphere, 
i.e., the angular correlation function, w(9,m), for an- 
gle 9 and survey "limiting apparent magnitude" m, 
can be more easily measured than £. This means that 
£ can be effectively measured at larger redshifts than 
is obtainable in redshift surveys — but at the cost of 
having to deduce £ from its sky projection. 

Hence, investigation of the z dependence of £ can 
be carried out to larger redshifts by measuring w 
rather than £. Many estimates of the amplitude of 
w(9, m) of faint galaxies in faint magnitude ("deep") 
small angle (a few square arcminutes, "pencil beam" , 
or up to a few square degrees) surveys have been car- 



small scale end (but noisy). 



' joveday et al. (1992) find their power law fit to £ 



200k- 1 kpc < r < 20/i _i Mpc 

Direct estimates of the evolution in this power law 
at low z me( j (m edian redshifts) include e = 1-6 ± 0-5 

0-4) 



et al. 1993; Roukema & Peterson 1994; Infante & 



ried out recently (e.g., Efstathiou et al. 1991; Ncu 
schaefer et~aT| 1991; jpritchet fc Infante] 19 92; [Couch] 



of t Ee Stromlo-APM rcdshift survey to be valid over Pritchct] 1995; |Braincrd et alj 1995). These angu 



Warren et al. 1993, 2 me( j= 
It al.| 1994 z med = 0-16) Q and e 



2-0±2-7 (Cole 

+ 1.0 







-1.3 



(IShep- 



icrd et alj 1996, z mec j= 0-36). However, consider- 
ation of galaxies as a single population, i.e., adop- 
tion of 7'o = Mpc and the median redshift 
cz mpA= 15, 200 km s -1 ) of the Stromlo-APM survey 



( Loveday et al. 1992, 1995) for comparison with the 
higher z mot \= 0-56 CF RS estimate of £ (Le Fevre et 
al. 199 6 ICFRS- VIII| ) would imply that e = 2-8 (cf. 
§4.1(3), pFRS VIII| ). 

This latter value of e is considerably higher than 
expected either for clustering fixed in comoving coor- 
dinates (e = 7 — 3 w —1-2, so that the index in Eq. [j]is 
zero); clustering fixed in proper coordinates on small 
scales ["stable clustering", e — 0, since the numbers 
of "clusters" changes by (1 -I- z) 3 and the number of 
galaxy pairs by (1 + z) 6 , the factor in Eq. |l| for r 
in proper coordinates is (1 + z)~ 3 ]; or linear growth 
of density perturbations in an Einstein-de Sitter uni- 
verse (e = 7 — 1 ~ 0-8). [Note also that the high 
z behaviour of £ may be different (e.g., Ogawa et 
1997).] 

A lower value of e could be justified by hypothesis- 
ing major changes in visible galaxy populations from 
low to high redshifts, as suggested by indications of 
colour dependence in the correlation function (e.g., 
infante fc Pritchct] 1995). 



lr The more recent estimate of Tucker et al. (1997) finds similar 
behaviour on scales down to about 20h~ 1 kpc, but for the 
redshift-space correlation function rather than the "real" space 
correlation function, which makes this difficult to interpret. 
2 Cole et al. 's (1994) e relates to e as e 



: 7 (l-e)-3. 



lar correlations are usually interpreted by integrating 
£(r, z) (Eq. |l|) over an appropriate domain in z — u 
space (see Eq. ||). 

The resulting values of e, or hypotheses about how 
population transitions might justify change in the 
value of ro as a function of z, are generally discussed, 
and an indication of the scale of galaxy pair sepa- 
rations is sometimes indicated by r±_ = dp IO p9 = 
r{z, u = 0) (for r in comoving coordinates; u is red- 
shift separation of a galaxy pair at mean redshift z, 
defined below for Eq. |2|) at the median redshift, z me( j , 
of the redshift distribution of the galaxy sample. 

However, the integral also includes galaxy pairs at 
unequal redshifts {u ^ 0) and galaxy pairs at mean 
redshifts lower and higher than z mec j, so rather than 
the usual assumption that r±_ is representative of typ- 
ical scales contributing to w, it is preferable to analyse 
the integral more closely. 

This is the goal of this article: to determine what 
pair separations contribute to w(9, m) over a region of 
parameter space judged likely for observational values 
of w(9, m). This is done by (a) separating out the dif- 
ferent contributions in the double integral relating £ 
to w, and by (b) defining an effective separation 
to be the median separation of galaxy pairs which 
contribute to a fraction p of the numerator of this in- 
tegral. The effective separation depends, in principle, 
on both the redshift distribution and the evolution 
of $, so r^g- is evaluated for likely ranges of relevant 
parameters. 

Indeed, the assumption that r± is a typical scale 
of galaxy pairs contributing to w turns out to be a 
good intuition, to better than an order of magnitude. 
In this paper we present quantitative justification for 



this intuition. 

In §p|, the double integration of £(r, z) is presented, 
using Sawicki et al.f s (1997) photometric redshift dis- 
tribution for the Hubble Deep Field (HDF; |Williams 
et Q 1996) as an illustration. The effective separa- 
tion, r^g, is defined and evaluated in In impli- 
cations of the resulting r^g- values are discussed and 
§|| presents the conclusions. All discussion is in co- 
moving (t = to) units unless otherwise noted and the 
Hubble constant is lOO/i kins" 1 Mpc -1 . 

2. LIMBER'S EQUATION 

The angular correlation function (small angle ap- 
proximation) is given by the the double integration of 

t{r,z), 



w(9, m) 



J dz N z (z,m) 2 J du£(r,z) 
[JdzN z (z,m)] 2 



(2) 



where 9 is the angle on the sky, m is the appar- 
ent magnitude, z — (z\ + Z2)/2 and u = z\ — Z2 
parametrise the redshifts of two galaxies at redshifts 
z\ and Z2, r{z,u) is the spatial separation between 
the two galaxies, and N z (z ,m) = d 2 N/d mdz is the 
redshift distribut ion at m (Limber 195 3; [Phillips et 
|l] 1978; |Pccbles| 1980; |Efstathiou et al.| 1991). 

One can think of £(r, z) as the excess probability 
that two randomly chosen galaxies lie at a (comov- 
ing) separation r. By symmetry, all such galaxies sep- 
arated in projection by 9 can be thought to form the 
surface of a cone of opening angle 9/2. Again by sym- 
metry, the dependence on the second angle </> can be 
dropped, so that integration is only needed over the 
double range of redshifts of pairs of galaxies, weighted 
by the numbers of galaxies at each of the two red- 
shifts. Note that the angular correlation function w 
is independent of the normalisation of d 2 N/ (dzdm), 
so only its shape is relevant to the integral (e.g. Yoshii 
1993). 

In order to understand the relative importance of 
the different factors in the integral, we separate them 
out in Fig. [I]. The redshift distribution chosen for 



this illustration is Sawicki et al.j's (1997) photometric 



redshift distrib ution for the HDF galaxies (see also 
Lanzetta et al. 1996) in the three Wide Field Camera 
images of the Version 2 (1996 February 29) drizzled 
version of the data to Iab,si4o < 27. The square of 
this distribution is shown in Fig. |](a). 

The redshift and distance dependence of £[r(z, u), z] 
can be seen in Fig. 0(b). At any given z, £ decreases 
as u increases, since the spatial correlation decreases 



for increasing r, quickly approaching zero. Both the 



use of a fixed angle 9 (chosen here as 9 = 4" 



since 
s w 



the range of significant signal in Villumsen et al 
HDF estimate is roughly 2-5" to 10") and the choice of 
"stable clustering in proper coordinates" (i.e., e = 0) 
contribute to the decrease in £ along the path for 
which u = 0. For large (fixed) it, r(z, u) becomes dom- 
inated by a nearly radial distance separation. Since 
d 2 dprop/dz 2 < 0, these separations become smaller 
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Fig. 1. — Analysis of how the spatial correlations £ com- 
bine with the number distribution dN/dz to form the 
integral w(9,m). Each of the panels is plotted against 
the same redshift range horizontally. From bottom to 
top, the four panels show: (a) [N z (z, m)Az] 2 in units 
of (squared) observed numbers of galaxies in bins of 0-2 
z units; (b) contours of log 10 {£[r(z,u), z]} in the z — u 
plane, showing the highest values at low z and low u; 
(c) contours of constant (dN/dz) 2 ^[r(z, u), z] in the z — u 
plane within which a fraction p of the numerator of Eq. ^ 
is obtained, for p — 10%, 30%, 70% (thin curves) and 
p = 50%, 90% (thick curves); (d) comoving separations (in 
kpc) corresponding to the u range plotted in (b),(c) — 
log w r(z,u = 0) and log 10 r(z,u = 4-lxl0 -4 ) are the 
lower and upper curves respectively. The redshif t distri- 
butio n is the HDF photometric z distribution of Sawick: 
et al. (1997). (The published histogram is smoothed by a 
gaussian of width a z = 0-1.) Parameter values are Qq — 1, 
A = 0, 9 = 4", ro = 5-5/2" 1 Mpc, 7 = 1-8, e = 0. 



for increasing z, which would imply higher £ for clus- 
tering fixed in comoving coordinates (e = 7 — 3) as 
z increases. Indeed, Fig. 0(b) shows that for "large" 
constant values of u, the increase in £ as z increases 
is strong enough to overcome the weakening of the 
amplitude of £ for higher z. 

Fig. ^(c) shows that neither one nor the other of 
these limits in u is sufficient to cover the domain in 
z — u space which is relevant for the total value of 
the double integral. The maximum values of u neces- 
sary to cover 90% of the numerator in Eq. are less 
than about 0-05% of a redshift "unit" , but this is still 
large enough that the dominance of radial separation 
over perpendicular separation causes £ to increase as 
z increases at a fixed it, at least for 9 — 4". 

The smallness of the value of u required to account 
for even 50% of the numerator of Eq. is the key ele- 
ment required to interpret observational estimates of 
w(6,m). In Fig. 0(d), the u range covered in the z — u 
plane plots is shown in units of (comoving) /i -1 kpc. 
Nearly all of the 90% contour is for galaxies separated 
by less than 1ft, -1 Mpc; much of the integral comes 
from separations less than IOO/2. -1 kpc. This is the 
scale on which galaxy interactions and the relation- 
ship between galaxy haloes dominated by dark matter 
and the visible galaxies containing stars are likely to 
be complicated. 

Other points to be noted in Fig. p](c) are: 
(1) 



Sawicki et al 



s (1997) redshift distribution 
includes many galaxies in the lowest redshift bins. 
Whether this is due to a genuinely steep faint end of 
the galaxy luminosity function or to the use of photo- 
metric rather than spectroscopic redshifts, it is not of 
cosmological interest to include spatial correlations of 
local dwarf galaxies within a few hundred kiloparsecs 
of the observer in estimates of w(z mec j) >, 1. It is there- 
fore of importance to observe that the percentile con- 
tours concentrate towards lower and lower z for higher 
values of the integrand (lower percentiles). The com- 
bination of rapidly decreasing separations, increas- 
ing values of £(r, z) and slowness of N z (z,m) to de- 
crease imply that a significant fraction of the signal 
in w(0,m) could be due to very local galaxies. 

For calculations presented below, a (conservative) 
lower limit is therefore set such that redshifts for 
which the perpendicular separation r(z, 0) is less than 
lO/i -1 kpc are excluded from the domain of integra- 
tion. 

(2) |Sawicki ct all 's (1997) N z peak at z « 2-3 con- 
tributes significantly between the 70% and 90% con- 
tours, i.e., roughly 10% of the numerator in Eq. p| is 



due to this peak. Given the uncertainties in Villum- 
sen et al.'s (1996) estimate of the HDF angular cor- 



relation function (essentially Poisson due to the small 
numbers of objects), this is not likely to be important 
for interpretation of the HDF data. However, for the 
brighter and more precise angular correlation func- 
tion estimates (e.g., Brainerd et alj 1995) this could 
be more important — analytical single-peaked redshift 
distributions may not be precise enough. 
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Fig. 2. — (a) (top panel) Dependence of r^g. [the median 
(comoving) radius which contributes a fraction p of the 
double integral determining w(6,m)] on the growth rate 
e of the spatial correlation function. The upper (lower) 
six cu rves are for p — 90%; (p = 50%). Curves are thick 



for sawicki et al.'s (1997) HDF 7as,si40 < 27 photometric 
redshift distribution and thin for the analytical distribu- 
tion of Eq. H with 20 = 1-20 (having the same z mec j as 
Sawicki et af['s redshift distribution) . Different cosmologi- 



cal models are shown [in (a)-(c)] by solid (flo = 1, Ao = 0), 
dotted (f2 = 0-1, Ao = 0) and dashed (fi = 0-1, A = 0-9) 
curves, (b) (middle panel) Dependence of r^g. on 20, 
adopting Eq. (for which 20 ~ 0-95 z me( j)) shown di- 
vided by minimum perpendicular distance at a given red- 
shift r± = r(2o,0). Upper (lower) three curves are for 
p — 90%; (p = 50%). (c) (bottom panel) Dependence of 
r± = r(zo,0) on 20 (equivalent to the angular diameter- 
redshift relation). Parameter values held constant in all 
curves are 8 = 4", ro = 5-5/i _1 Mpc and 7 = 1-8. 



3. THE EFFECTIVE SEPARATION r p „ 

eft 

Of course, the values of the separation scales which 
contribute most of the value of w(9,m) discussed 
above are for a particular choice of the parameters 
chosen in Eq. |l| (or on the validity of this equation 
at such scales), the redshift distribution, 9 and the 
cosmological model. In order to discuss dependence 
of the separation scales on Eq. []], on the z med and the 
shape of the redshift distribution, on and on cosmo- 
logical parameters, it is useful to define the effective 
separation, r^g, to be the median value of r(z, u) over 
the domain 'm z — u space inside a contour of constant 
N z (z, m)t;[r(z, u), z] which contributes a fraction p of 
the value of w(8, m) (where the denominator in Eq. ^ 
is held fixed). 

3.1. Dependence on the Spatial Correlation 
Function 

Lowering the value of ro (which is considered by, 
e.g., Bernstein et al. 1994; Brainerd et al. 1995, to 
denote a possible transition in galaxy "populations" ) 
factors out of Eq. |as a constant; it reduces the value 
of w but does not affect the question of which separa- 
tions are relevant to the integral. Reducing the value 
of 7 (for fixed e — 7 in order not to affect the red- 
shift evolution of £) would increase the values of r^g 
slightly. In estimates of w for faint galaxies, it is usu- 
ally assumed that 7 — 1 = 0-8 in order to correct for 
the integral constraint, so this is the value of interest 
here (e.g., Hudon fc Lilly 1996; though note Campos 



it al. 



1995; [Infante fc Pritchet| 1995). 
Of more interest is the value of e, which as men- 
tioned above has either theoretically or observation- 
ally motivated values lying in the range from —1-2 to 
nearly 3. 

Fig. ||(a) shows that in spite of the wide range in 
values of e, the dependence of r^g on e is weak. Even 
for a cosmological constant do minated flat co smolog- 
ical model and e = 3, r^g for 3awicki et al.'s (1997) 
redshift distribution is no more than a factor of two 
lower than that for the analytical distribution. 

Although weak, there is an average trend to lower 
values for higher e. Higher e means a more rapid de- 
crease in spatial correlation amplitude with z, so that 
the correlations at lower z contribute relatively more 
to the integral. The lower the value of z, the lower 
the (comoving) perpendicular separations r(z,0), for 
any cosmological model. Hence, the slight decrease in 

r eff 

In other words, independently of the growth rate of 
the spatial correlation function, angular correlations 
in deep surveys of faint galaxies only probe length 
scales characteristic of galactic haloes. 

The stronger decrease in r^g for the observational 
(photometric) redshift distribution than for the an- 
alytical distribution of the same z med can be at- 
tributed to its complex shape. Once e has increased 
enough that the z ~ 2-2 peak no longer contributes 
much, the effective z med of the observational distri- 



bution will be lower than that of the analytical dis- 
tribution; thus, the lower r^g values. 

The trend to lower r^g is not monotonic for the ob- 
servational redshift distribution. This is attributable 
to an increase in the relative importance of u separa- 
tions to perpendicular separations when N 2 (z, m)(l+ 
z )-(3+ £ t) becomes flatter. A "complex" enough red- 
shift distribution could possibly increase the impor- 
tance of this reversed dependence of r^g on e, as in- 
deed is seen below. 

3.2. Dependence on the Redshift Distribu- 
tion 

The effective separation r^g depends on a redshift 
distribution primarily via the typical rcdshifts in the 
sample, or more quantitatively, on the median red- 
shift, 2 mec i, of N z . This dependence is analysed here 
by use of the simple analytical distribution 



z 2 exp[ 



-0/%) 



(3) 



where (3 — 2-5 (as in Efstathiou et al. 1991; Villum- 
5cn et al. 1996). For this distribution, z mec j depen- 



dence translates directly into dependence on zo, since 
z med w 0.95z - 

A very useful property of this dependence is shown 
in Fig. ||(b): nearly all of this zq dependence is in the 
(equivalent of the) angular diameter-redshift relation, 
shown as the dependence of the perpendicular sepa- 
ration r±(z) = r(z, u — 0) on z in Fig. ||(c). 

Both r^g" and are nearly constant over a 

wide range in zq. Moreover, r 5 g % is only about 5% — 
15% greater than r±(z ), i.e., about 10% — 20% 
greater than r±(z med ). Relative to the precision of 
present cosmological measurements, it can be said 
that r±(z mec f) is a very good estimator for r^° for 
such analytical distributions. In addition, since r^° 
is insensitive to e and since the e dependence of both 
the observational (photometric) and analytical N z 
curves are similar (Fig. |^), the quality of ^i(-z me d) as 
an estimator for r^g° is little dependent on the value 
of e and the precise shape of N z , avoiding the need 
for explicit integration of Limber's Equation (Eq. ^|). 
The 90% percentile median separation 

r eff % is 

about 4 — 4-5 times r±(zo), depending on zq and cos- 
mological model to no more than about 10%. Again, 
the relative robustness of r^g to e implies that a con- 
stant value of ^gg % in this range can be used for prac- 
tical purposes without integration of Limber's Equa- 
tion being required. 

The dependence of r P g on the shape of N z is not 
totally negligible. This can be seen to some extent 
in Fig. ||(a) , which shows that for the the 
values of r^g can differ by up to about a factor of two. 

To confirm this more explicitly, various published 
photometric redshift distributions estimated for the 



HDF have been adopted in place of that of Sawicki 



ct al. (1997) and the resulting dependence of r P g on 
e for these N z estimates are presented in Fig. ^ The 
N z estimates adopted are shown in Fig. [|. 




The results are similar to those for 3awicki et 



Fig. 3. — The dependence of r^g on the shape of jV ; 
is shown by reproducing Fig. ^](a), but repla cing the N, 
of ISa wicki etafl (1997 ) by those of (a) 
Lib " 



Lanzetta et 



estimate of 

al.[ ]l996), ( b) |Mobasher et al.| (1996, Fig.2) and (c) |Gwyn 



fcl |artwick| (1996), and using analytical distributions with 
appropriate 20 values (a) 20 = 1-15, (b) 20 = 1-6 and (c) 
2 = 2-0. Axes and line styles are as for Fig. H(a). 



-1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — |- 
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Fig. 4. — The shapes of the photometr ic redsh ift distri- 
butions N z esti mated by Sawicki et al. (1997), [Lanzetta 
et p] (1 996), |Mobasher et al.| (1996) and pwyn &: 



Haitwick (1996) for the HDF are shown here against red- 
shift, normalised for purposes of comparison. The depen- 
dence of r^g. on the different shapes is shown in Figs 2(a), 
3(a), 3(b), 3(c) as indicated. 



s (1997) N z , though the maximum difference be- 
tween photometric and analytical distributions in- 
creases from about a factor of two to a factor of three 



to four for 



„90% 



modality of Gwyn fc Hartwick 



in the case of h igh e values 
s 



high resolution of Mobasher et al 



The bi- 
(1996) N z and the 
(1996) N z make 



this hardly very surprising. 

The large differences between the photometric and 
analytical distributions for Gwyn fc Hartwick 's (1996) 
N z are mainly due to the low z peak. Obviously 
there is uncertainty in the correctness of photometric 
redshifts, but if N z (HDF) really was as bimodal as 



Gwyn fc Hartwick (1996) estimated with a strong low 
z peak, then not only would many galaxy pairs which 
contribute to wo be closer to one another than would 
be expected from an analytical unimodal distribution, 
but a large part of the wo signal would be from non- 
cosmological distances, at which r± ~ IO/1 -1 kpc. In 
this case, the meaning of a wq estimate would have 
to be analysed in considerable detail, or else redshift 
separation would have to be used to estimate w(r p ) 
(Davis & Peebles 1983) rather than wo in order to 
extract quantities with simple physical meaning. 

In summary, the earlier photometric redshift dis- 
tributions to that deduced by Sawicki et al. (1997) 
imply similar values of r^g to those of Sawicki et al 



to within much less than an order of magnitude. 

It should also be noted that changing the limiting 
magnitude of the survey or changing the value of (3 is 
equivalent to a change in z mec n and/or in the shape 
ofN z . 

3.3. Dependence on Other Parameters 

Since r{z, u) is dominated by the perpendicular 



(z) (for 9 < 1 rad), 



separation r{z,u = Q) = 9 dprop 
r^g. is nearly proportional to 9. Indeed, although the 
line-of-sight galaxy separation for a given z is not di- 
rectly proportional to 9 [i.e., r(z,u,92)/r(z,u,9i) ^ 
O2/O1 if u 7^ 0], the required domain in u to obtain 
a fraction p of the integral of w(9, m) expands suffi- 
ciently (but remains small enough) that the value of 
r^g- remains proportional to 9 to the precision of the 
calculation (for p — 50% and p — 90%) . 



Values of 



„50% 



Villumsen et al 



r 90% 
eff 



for the nominal angle of 



s (1996) w determination [w(9,m) 
is normally interpolated or extrapolated to the same 
angle for a large number of different surveys in order 
to compare the relative amplitudes of the surveys], 
i.e., 9 — 1", imply r V rr values four times lower than 



those displayed in Fig. 0, i.e. 



,90% 
r eff 



r 50% 
r eff 



lO/i 1 kpc and 

40/i _1 kpc. While such a nominal angle is 
primarily meant for comparison of amplitudes of w, 
the physical meaningfulness or otherwise should ob- 
viously be kept in mind when examining such figures! 

Brighter faint galaxy surveys having most of their 
signal at, say, 40", would have their r p „ values mul- 



tiplied by ten, 
1500/t- 1 kpc. 



r 50% 
eff 



eff 

SQOh" 1 kpc and r 



90% 

eff 



The effects of cosmological parameters on r P g- are 
shown in panels (a) and (c) of Fig. || since ddp m p/dz 
is successively larger for hyperbolic and cosmologi- 
cal constant dominated, flat cosmologies, the values 
of r{zo, 0) and consequently r P g. increase respectively. 
The increase in r(zo,0) is in fact faster than that of 
r P g, so that Fig. |](b) shows a slight reversal in this de- 
pc ndcncc for the ratio r p fF /r(zp, 0), though the main 



effect is in t(zq, 0). 



4. DISCUSSION 

The results above can be summarised as order of 
magnitude estimates: 



v 



h3ff- 



(1 +A p )r ± [* mcd (/V a )] (4- 



where 



r \ z = z med( N *)> U = 0, fi = I, 
A = 0] 

fldproptzmedC-Nz)'^ = 1,A = 0] 
2c „ A 1 \ 



Ho 



1 



and 



1< A 50 <2 
2 < A 90 < 7 



(5) 



(6) 



for 7 = -1-8, -1-2 < e < 3, 0-1 < O < 1-0, 0-0 < 
A < 1 — f2oi 1 < ^medl-^z) < 2-5 and for N z shapes 
as "smooth" as in Eq. [| or no "rougher" than that of 
sawicki et al.'s (1997) redshift distribution. 



These relations enable r p a to be estimated for 

eft. 

observational angular correlation function measure- 
ments. Many of the brighter (B < 24, V £ 25 or 
R < 26) faint galaxy angular correlation measure- 
ments (|Ncuschaefcr ct al. 1991; Pritchct & Infante 1992; 



Couch et al. 1993; Infante fc Pritchet| 1995) show sig- 
nificant power law correlations on scales ranging from 
~ 10" to ~ 10'. Adopting z me{ j ~ 0-5 as a typical 
value for the above surveys, which should give r± 
to be correct within better than an order of mag- 



nit 
lat 



do, Eqs || (which remain valid under extrapo 



ibs- 



to z 



mtsd 



0-5) imply that bQh^ 1 kpc 



50% 



,3/j- 1 Mpc and 200/1" 1 kpc 



„90% < 
r eff 



efr 



12/1" 1 Mpc 



for these angular ranges. Spatial correlations are well 
established at most of these spatial separations, so in- 
terpretation in terms of Eq. [j] is likely to be a good 
first approximation. 

Possibly of interest in these cases is the higher end 
of the length scales. Loveday et al.| (1992) only con- 
sider their power law fit to £(r) to be good to about 
20/i~ 1 kpc. The surveys covering the largest solid an- 
gles may well be affected by — or be used to estimate — 
spatial correlations at separations difficult to measure 
by local surveys such as the Stromlo-APM survey of 



Lo veday ct al. Indeed, Infante & Pritchct (1995) di 



large scale behaviour of £(r). Their angular correla- 
tion functions extend to about 0-4°. The analysis here 
implies that at 0-4°, ~ 30/1" 1 Mpc, so careful 

analysis of this data might be able to produce esti- 
mates of the perturbation power spectrum on scales 
comparable to that in the Stromlo-APM survey. 

The brighter ang ular correlation m easur ements on 
smaller field s (e.g., Efstathiou et al.| 1991; Roukcma 



& Peterson| 1994; |Brainerd et aL 1995) cover smaller 
angular ranges, from ~ 10" to ~ 1', so again adopting 

,90% 



,.50% 

r eff 



< r 



off 



z med ~ 0-5 implies that 50/i 1 kpc 
<,1.2/i _1 Mpc for these observations. This is again 
within typical estimates of power law behaviour for 
£(r), e.g., Davis fc Peebles| (1983) and Loveday et 
al. (1992), though 50ft. -1 kpc is already smaller than 
typical estimates of the radii of the haloes of L* 
(bright) galaxies. 

Raising the es timate of z mp c \ from 0.5 to z me( ^ = 
1-2 (the value for [Bawicki et al. 's HDF analysis) would 
double each of these separation sizes. Alternatively, 
conversion to proper units, which are relevant if con- 
sidering galaxy haloes at z « 1 to be collapsed 
and dynamically stable objects, requires division by 
(1 + z) ~ 2, i.e., halving the separation estimates. 

This brings us back to the worry raised above — 
that spatial correlations are being measured for galax- 
ies which are closer together than the radii of their 
parent haloes! However, since a large part of the 
signal in the cases just-mentioned comes from larger 
separations, this should not affect the observational 
validity of using Eq. | as a first order estimator. 

The reliance upon small separations becomes a 
stronger concern for the angular correlation function 
of the HDF galaxies (and for possible future estimates 
of w(0, m) using the HST, the Next Generation Space 
Telescope, or if adaptive optics can be used for faint 
fields using ground-based telescopes) . 

The I A b < 27 limit of |5awicki et alj 's (1996) anal- 
ysis of the HDF galaxies corresponds to 26-5 (in 
the standard Vega-based system), and t he range of 
the galaxies' colours (e.g., Metcalfe et al.| 1996) is ap- 
proximately $ R — 1^2. He nce, the redshifts ap- 
propriate for Villumscn ct al.'s (1996) i?-band mag- 
nitude ranges R < 26-5 to R < 28-5 used in estimat- 
ing w(9,m) should be similar to those in Bawicki et 



al. 



s analysis, though the bluest galaxies not present 



in Sawicki et al.'s (1997) analysis could be expected 
to be at somewhat lower redshifts. If we consider 



the range of significant signal in Villumscn et al 



(1996) estimate to be from 2-5" to 10", this then im- 



,,50% 
r eff 



< r 



.90% 
eff 



plies for z mec j « 1-2 that 25h~ 1 kpc < 
<, 400/i -1 kpc for these observations. 

In this case, almost none of the signal comes from 
separations for which Loveday et al. (1992) found sig- 
nificant correlations fit by a power law, and the sim- 
ilarity to halo sizes is significant. Possibly one of the 
neatest measurements of halo extent is by Lyman-a 
and metal line absorption systems in front of quasars. 
For instance, Bergeron fc Boisse] (1991), Bechtold ct 



cus b the degree to which their data io affected by the 



al. (1994), Lanzetta et al. (1995), Fang et al. (1996) 



and Le Brun et al. (1996) estimate gaseous halo radii 
for both sorts of systems of around 50 — 200ft- 1 kpc. 
This range matches very closely the range in r^g just 
mentioned. 

[However, it should be noted that (1) the Lyman- 
a systems may not be associated with galaxies (e.g., 



However, Davis & Peebles's (1983) observation of 



Le Brun et al. 1996) and (2) there is marginal rota- 
tion curve evidence that the h alo of the Galaxy is no 
more t han 15 kpc in rad ius ( Honma fc Sofue 1996, 
but sec Binncy & Dchncn 1997), implying a mass for 



the Galaxy about an order of magnitude smaller than 
most estimates.] 

Is the value of z mec j adopted here too low? While 
gawicki et alj 's (1997) photometric redshift analysis 
uses template spectra which account for both inter- 
nal reddening and high-z Lyman absorption, earlier 
analyses, e.g., Mobasher et al. (1996), suggested that 
z med should be as high as £ me( j = 2-1. Even if the 



red shifts of th e i?-selected sample contributing to Vil 
lun .sen et al. 's w(9,m) estimates (or to independent 
estimates of, say, an /-selected sample) were as high 
as z rs 2-1, these effective separations would only 
increase by about 30% — 50% (depending on cosmol- 
ogy), according to Eqs H-pl The spatial separations 



would still correspond to typical estimates of galaxy 



halo radii, particularly when applying the conversion 
from comoving to proper coordinates. Note that the 
values adopted for use in Eq. [| in Villumscn et al.'s 
(1996) models shown in their Fig. 2 for R < 26-5 to 
R < 28-5 are 1-4 < z < 1-8. 

This result brings to mind the suggestions, based 
on morphological analysis of the HST Medium Deep 
Survey (MPS) a nd th e HDF QGriffitlis et aE| 1994 ; 
Cascrtano et al.| 1995; Driver et al. 1995a, b; Glaze- 



brook et al. 1995; Abraham et al. 1996; van den Bergh 



et id] 1996), that many of the HDF galaxies are the 
"building blocks" of future galaxies which later merge 
together. ( Pascarelle et alj (1996) make a similar sug- 
gestion based on HST imaging and Multiple Mirror 
Telescope spectroscopic redshifts.) In this case, the 
HDF galaxies' haloes should be smaller than present 
day galaxy haloes, and so the visible galaxies could 
be correlated without their haloes necessarily overlap- 
ping. 



Table 1: Values of r f ' a for 1" for redshifts z = 2-5 and 
eft 

z = 5 spanning redshifts likely for the HDF objects 



col our-selected to high redshift by |Colley et al. (1996 



19S 7), assuming that Eqs 4-3 can be validly extrapo- 



lated to these objects. These effective separations are 
in units of ft -1 kpc, comoving or proper as indicated. 



comoving 



proper 





n 


Ac 


z = 2-5 


z = 5 


z = 2-5 


z = 5 


50% 


1-0 


0-0 


14 


17 


4 


3 


50% 


0-1 


0-0 


21 


35 


6 


6 


50% 


0-1 


0-9 


26 


37 


7 


6 


90% 


1-0 


0-0 


54 


69 


16 


12 


90% 


01 


0-0 


86 


141 


24 


24 


90% 


0-1 


0-9 


103 


147 


30 


25 



correlated galaxies down to lOft -1 kpc is a local obser- 
vation, it is not an observation of primordial galaxies. 
It is quite noisy, so perhaps is already affected by the 
detailed nature of the galaxy-halo relationship. In ei- 
ther case, the angular correlation function of the HDF 
must be affected by, and hence offer some clues to, the 
relationship between galaxies and haloes. 

A compl ementary ana lysis to that pres ented here 
is that of ICollcy et al] (1996; 1997). |Colley et 
al. (1996) measured a very strong angular correlation 
for "galaxies" selected by colour to be at high red- 
shifts [z > 2) and point out that the detected objects 
may in fact be star-forming regions within "normal" 
galaxies, dimmed only by (1+z) 2 in surface brightness 
rather than (1 + z) 4 since they are (nearly) unresolved 
by the HST. In |Colley et al"] (1997), apart from a dy- 
namical discussion of different scenarios, the authors 
add a nearest-neighbour analysis which strengthens 
their argument that there is an excess of close objects 
separated by less than about 1". 

The result of the effective separation calculations 



presented here only strengthens Colley et al. ! s model. 
As they point out in their conclusions flColley et 



ah 1997), two galaxies at identical redshifts sepa- 
rated by 1" at the redshift range of their sample are 
separated by roughly 6 kpc in proper units. Since 
r^g- ~ rj_, if the objects were indeed individual galax- 
ies, and if a much larger field were measured in order 
to reliably measure w(9, m) at these small angles and 
faint magnitudes, then about 50% of the spatial cor- 
relation signal would be from galaxies separated by 
less than about 6ft -1 kpc, and 90% of the signal for 
galaxies closer to one another than about 24ft _1 kpc. 
These values, in both comoving and proper coordi- 
nates in order to facilitate deductions from both cos- 
mological and galaxy formation points of view, are 
presented in Table 0. 

Returning to the full sample of the HDF galaxies, 



for which 25ft. kpc 



„50% 
r eff 



< r 



,90% 

eff 



,400ft- 1 kpc, 



are there any justified alternatives to the extension of 
Eq. |l| to these separations, either observationally or 
theoretically? 

A comparison of the theo retical evolution of the 

's (1996) 



spatial correlation function to Villumscn et al 



HDF estimates is given by Matarrese et al 
for an Qq = 1,Aq 



(1996), 



= CDM universe where the 
full linear and non-linear evolution is given by a for- 



mula fit to gravity-only N-body si mulations ( Hamil 

~ 1995) 



ton et al. 1991; Jain et al 



Matarrese et al 



(1996) Fig. 5 (for no bias factor) is approximately 



consistent with the HDF correlations, whilst Villum 



sen et al. (1996) found that a spatial correlation 
growth rate of e = 0-8 in Eq. ^ was not sufficient 
to provide low enough correlations unless the value 
of ro is decreased, i.e., that a transition in popu- 
lations is hypothesised. Differences at the sub-Mpc 
level could explain this better fit. However, the CDM 



simulations used for Jain et al.'s formula have only 
1 particle/(350ft _1 kpc) a and a force resolution of 
32 ft -1 kpc, so may not be valid at these scales. Ad- 



ditionally, of course, they cannot address the ques- 
tion of how to relate stellar galaxies to dark matter 
dominated haloes — unless this is governed by gravity 
alone, which seems unlikely. 

An empirical, but theoretically motivated, alter- 
native would be Peacockj 's (1996) analysis based on a 
two-power law fit to galaxy correlations of the APM 
( |Maddox et"aT] 1990a; 1990 b) and IR AS ([Saunders et 
al.|Jl992) surveys. However, [Peacock| finds that this fit 
implies that for an CIq = 1, Ao = universe, an extra 
correlation at scales of r < 2irh~ 1 Mpc is needed to 



fit the CFRS flCFRS-VlH 1996) data, unless a scale- 
dependent bias factor is used in combination with a 
constant bias factor. Matarrese et al.'s (1996) use 



of Jain et al. 's (1995) N-body motivated £(r,z) de- 
scription better fits the CFRS data (for no bias), so 
would maybe more appropriate for an extension of 
the present work. 

An alternative empirical improvement to Eq. |l| 
could be to use the w estimates of Infante et aL (1996) 
on very small scales (down to w 1") in a survey cover- 
ing a very large solid angle (2-3 sq.deg.), hen ce, hav- 
ing enough precision to make such estimates. Infante 
et ^l] (1996) find that at the smallest separations, 
2" <6 <6", there is a correlation signal about a fac- 
tor of three higher than that of a j = 1-8 power law 
extrapolated from larger angles. To the extent that 
the above calculations can be extrapolated to this 
"highly non- linear" case, the value r± — 24ft, -1 kpc 
(as defined in Eq. |J) implies that half of this signal 
comes from three-dimensional galaxy pair separations 
of about this size, and 40% more comes from sepa- 
rations up to about 100 — 200/i _1 kpc. (Note that 
infante et al. adopt proper units. The conversion fol- 
lows from using z mec j = 0-35.) The values of £(r, z) for 
r < 20h~ 1 kpc, without contamination from "normal" 
correlations at larger angles, could therefore be even 
higher than might be expected from simple inspection 
of Infante et al.'s (1996) figures. 



5. CONCLUSIONS 

In order to see what length scales and redshifts 
correspond to the spatial correlations £(r, z) repre- 
sented in recent measurements of the faint galaxy 
angular correlation function w(6,m), Limber's Equa- 
tion (Eq. ||) has been examined over a range of likely 
possibilities for £ and N z (z,m). The separation into 
N z (z,m) and £[r(z,u),z] and the contours in the 
z — u plane which contribute to the integral have 
been shown, revealing that observational w measure- 
ments combine spatial correlations from a wide range 
in mean redshift. 

A simple definition of an effective separation r^g to 
quantify the relevant length scales for w has been de- 
fined as the median (comoving) separation within the 
contour of constant {dN / dz) 2 £ > [r{z, u), z] contributing 
a fraction p of the numerator of Eq. |[ This is similar 
in order of magnitude to the perpendicular distance 
r± = dpropO at the characteristic redshift, Zq, of sim- 
ple analytical redshift distributions (Eq. 0) , and hence 



to the median redshifts of such distributions. 

More precisely, for the parameter space covered, 
is about 5% — 15% greater than r±(z ), i.e., 



„50% 
r eff 



about 10% — 20% greater than r±(z me< ^) (for (3 = 2-5 
in Eq. ||); while is about 4 — 4-5 times r±(zo). 
For a more realistic redshift distribution, e.g., the 



photometric redshift distribution of Bawicki et al 
(1997), r^g. may be different by up to about a fac- 
tor of two to the r^g- values for the analytical redshift 
distributions if —1-2 < e < 3. 

In summary, values of r^g for typical median red- 
shifts of most faint galaxy angular correlation mea- 
surements indicate that the use of Eq. |l| does not im- 
ply an extrapolation (in galaxy pair separation) much 
below or above observational spatial correlation mea- 
surements. (The validity of redshift dependence scal- 
ing by (1 + z)~( 3 + £ ~t) i s no t tested, however.) 

The scales on which HDF galaxies are correlated 



(as measured by Villumscn et a\\ 1996) imply co- 



moving separations of 25h~ r ~ kpc r eff ° ^ r eS 
< 400/i _1 kpc. That is, for typical redshifts over z £ 1, 
a substantial fraction of the angular correlation signal 
is generated by galaxies spatially separated, in proper 
units, by around 10— lOO/i -1 kpc. While spatial cor- 
relations of galaxies at these separations have been 
observed for local galaxies (Davis & Peebles 1983), 
this is a scale at which halo and/or galaxy existence, 
let alone interactions, might strongly modify the spa- 
tial correlation function. 
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of the late Prof. Roger Tayler, whose support for 
this collaboration was well appreciated. We thank 
Alain Blanchard for useful discussion on this topic 
and B. Mobasher for providing tables in electronic 
format. BFR acknowledges a Centre of Excellence 
Visiting Fellowship (National Astronomical Observa- 
tory of Japan, Mitaka), support from PPARC and 
from Starlink computing resources. 
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